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I. INTRODUCTION 



There have been remarkable developments in our understanding of high-energy description of Quantum Chromo- 
dynamics (QCD) in the last decade. Although, in principle, main constituents of nuclei could be thoroughly governed 
by QCD dynamics, it would require enormous numerical computation to perform a first-principle calculation of QCD 
for the properties of nuclei. Most of problems encountered in QCD are non-perturbative, and besides, non-linear, 
meaning that gluons strongly interact with themselves. 

Under some conditions feasible in experiments, however, one may hope to derive analytical information directly 
from QCD in a manageable way, which is not only of general interest in theoretical physics, but also of importance 
as an experimental testing ground for theoretical ideas in QCD. It is well known that the strong coupling constant, 
as, is expected to be small for far energetic processes owing to asymptotic freedom, so that the perturbation theory 
would work. In fact, however, the perturbative expansion becomes reliable for processes in which the exchanged gluon 
carries a large momentum and it is possible that the momentum transfer remains small even for energetic processes. 
Moreover, even if the exchanged gluon conveys a sufficiently large momentum transfer, small as is not necessarily 
adequate to guarantee the validity of perturbative calculations. 

Even when the perturbative expansion in terms of as seems plausible, resummation is needed whenever enhancing 
logarithms exist in the expansion series, which occurs in general for processes characterized by multiple (large) energy 
scales. This is actually the case for the processes we are discussing in this paper in which partons with small Bjorken-x 
are relevant. Then, two large scales Aqcd ^ i <C s may give rise to a large logarithmic factor ~ ln[s/i] ^ ln[l/x]. 
Here Bjorken-x is intuitively interpreted as the longitudinal momentum fraction x = of a parton inside a 

target seen in the infinite momentum frame, — > oo; the target moves nearly at the speed of light in the positive 
z direction, that is, along the x+ — (z + 1)/^/2 direction in the light-cone coordinate. Even though as is far smaller 
than unity, the perturbation theory in as breaks down if aslii[l/x] turns out not to be small for x <C 1. In the 
leading-log approximation the resummation program with respect to such logarithmic factors has been developed, 
which brings about running in x for gluonic properties such as the gluon distribution. This is analogous to the 
famous (momentum transfer) running of gauge coupling constant after the summation over leading ln[(5^/AQQj-,] 
terms is taken in perturbative QCD. Insofar as only the perturbatively calculable part is concerned, it has been 
our common knowledge that such running in x, or the small-x evolution, in the leading-log approximation obeys the 
Balitsky-Fadin-Kuraev-Lipatov (BFKL) equation 0,0,13, H • 

The BFKL equation lacks the saturation effects |a,|a,l3, however, that would be important physics when the parton 
(gluon) distribution is highly dense in the small-x region. Since the BFKL dynamics encompasses the evolution in 
X but not in Q^, the transverse area is fixed and the number of partons increases in a fixed area as x goes smaller. 
Partons eventually overlap each other, indicating that non-linearity of QCD should be properly taken into account 
at a certain saturation scale Q^{x), and the parton number density must be saturated via recombination processes. 
The modern picture of this saturation mechanism is more suitably described by coherent gluon fields A'^ representing 
small-x wee partons 0, S 0- The larger-x partons are integrated out to lead to the gluon source p"^, from which 
the coherent fields are created corresponding to smaller-x degrees of freedom. [This is a similar concept to the 
Weizsacker- Williams approximation.] When A'^ becomes as large as ~ l/g, non-linearity is crucial for the saturation. 
In the first approximation the source is, originally in literatures, just regarded as the color charge density carried by 
valence quarks inside the target hadron, supposing that the target is a large nucleus. Since quarks in different nucleons 
are uncorrelated, the color charge density can be considered as distributing randomly. The random distribution of 
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strong gluon fields A'^ is commonly referred to as the color glass condensate [loj- 

One can compute the small-x evolution, that is actually the renormalization group equation with the rapidity 
variable, from x down to x — dx in the presence of the color glass condensate by assuming that one has such an 
effective theory of the color glass condensate at a certain x and one can retain the same framework at x — c?x after 
integrating out any degrees of freedom within the slice c?x = dp+/P+. The equation obtained in this way is to be 
considered as a natural extension of the BFKL equation with the saturation effect of dense gluons incorporated. 
There is an intrinsic energy scale (5g(x) which is large, so that the perturbative resunimation becomes infraredly 
stable The theoretical framework has been well organized and the central result is now widely known as the 

Jahhan-Marian-Iancu-McLerran-Weigert-Leonidov-Kovner (JIMWLK) equation |TollTllll[Tllll[Ta.lT7lfT^. 

The purpose of this paper is to present the derivation of the JIMWLK equation in a way different from the 
conventional arithmetic steps. The reasons why we shall attempt to rederive the well-established equation here 
are twofold: The author has recently proposed a representation for t he gaug e invariant source terms [l9j which is 
equivalent with the form as used in literatures conventionally [T^ Il4 lisl ITgL IitL ll^. but is rather suitable for 
treating the gauge invariance and non-commutative nature among color charge operators transparently. It would 
be a first-step check to confirm explicitly that such a representation proposed by the author really leads to the 
correct JIMWLK equation in the dense regime, which is necessary as a preparation for further applications. Another 
issue, that could provide a hint to approach the dilute regime where non-comniutativity among color charge density 
operators is important, is that in our formalism we can clearly see under what approximation we have disregarded 
non-commutativity in the dense regime. The JIMWLK problem, i.e., the derivation of the evolution equation in the 
dense regime, has already been well understood itself and our final results shall not add any new materials to that. 
Nevertheless, looking at the derivation of the JIMWLK equation from another point of view as addressed in this 
paper should be useful to understand the JIMWLK dynamics more deeply and would give us a clue on the treatment 
beyond the approximations utilized in the dense regime. In short, the ideology of this paper is; the derivation is new, 
while the derived equation is already known. 

In Sec. |n] we will quickly look over the philosophy of the color glass condensate and will articulate our final goal 
in advance, namely, the JIMWLK equation. We will present the central part of this paper in Sec. IIIII where we 
will explain a step-by-step derivation of the JIMWLK equation along the formalism involving quantum color charge 
density. Our summary and conclusions arc in Sec. IIVI Some discussions are followed in App. ^ to account for 
connections between our derivation and the conventional one for the sake of more clarity. 

II. COLOR GLASS CONDENSATE 

We shall reiterate the physical situation we are considering here, though we have partially explained it. A nucleus 
(target) is moving fast toward the positive z direction. In the light-cone coordinate the target is traveling along the 
temporal ~ {t + z)l\/2 direction while sitting on — if. — z)/y/2 ~ 0. In the momentum space language it bears 
a large longitudinal momentum p+ — {po +Pz)/ ^ oo. In the formulation of the color glass condensate the essential 
assumption is that there is a separation scale A for p+ ; all the fast partons with |p+ 1 > A are integrated out to result 
in the classical color source p%{x) which does not depend on x~^ (i.e. static) because of the time dilatation relative to 
slower partons. The color orientation is assumed to distribute randomly at each spatial point x = (x^,x±) with a 
weight yVA[/OA]- Under these assumptions the QCD generating functional is modeled as 

Z = JvpaWaIpa] J VA5[A+]exp{iSYM[A]+iSw[A-,pA]}. (1) 

It has been shown in that the distribution p%{x) is discrete at each spatial point at the microscopic level. In this 
paper, however, we shall work only in the dense regime so that we can approximate p^ix) as a continuous function 
as treated conventionally. Here we have employed the light-cone gauge (A+ = 0) that enables us to cut-off the 
longitudinal momentum in a way free from a mixture of momenta carried by A~^ through the covariant derivative 
(9+ — igA^. It is quite important to note that p%{x) has a support in the longitudinal direction as < < 1/A 
as a result of the integration over fa st g luons with an appropriate boundary condition imposed. This longitudinal 
structure is first recognized in |0, 0, 0| and can be also embedded in the weight functional Wa [pa] in a way leading 
to p%{x) — for a;^ > 1/A The JIMWLK equation is, as we will see later, consistent with the support structure 
as it should be; along with the quantum evolution in yVA[pA] the edge of the support diffuses to larger x~ . In our 
notation we do not have to attach the scale to the source specifically once Wa[pa] accommodates the support structure 
, but we will stick to denoting p\{x) in this paper to remind us of the longitudinal structure. 

The support in is intuitively to be interpreted as a coincidence of the momentum rapidity and the space-time 
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rapidity of fast partons as argued in That is, the rapidity variable, 



T = In 



(2) 



is defined as a momentum rapidity difference between the target and the parton. If the particle moves along the 
longitudinal direction and its velocity is not changed, which is actually the case in the eikonal approximation, the 
momentum rapidity is identical to the space-time rapidity. Therefore, if the target is at the position of Xq , then 
the parton with the rapidity difference r is located at x~ = x^e^ . Since p\{x) has an origin from fast gluons with 
T < ln[P"'"/A] in the color glass condensate picture, its distribution is within Xq ~ < x~ < 1/A ~ x~ . Hereafter 
we will often make use of r to mean a reference rapidity r — ln[P+/A] instead of using the separation scale A like 
Wri/Or] signifying >Va[pa] equivalently. 

In this paper we expressly call p\{x) the classical color charge density to distinguish it from the quantum one that 
fulfills the covariant conservation and the canonical commutation relations at the operator level. [We will explain 
what this means exactly in the next section.] In short p\{x) turns out to be not a source for quantum fields but 
a distribution of fast partons, which is a well-defined concept at the classical level. The reco gnit ion of p%{x) as a 
classical quantity characterizing the particle distribution has been discussed already in Ha 113, 13 where our p%{x) 
is denoted by i' in order to distinguish it from the genuine source for quantum fluctuations. Its significance has also 
been emphasized recently in |l9l |. 

The source action Sw , Pa] in is required to be gauge invariant and to be close to the simple eikonal coupling 
~ tr{pA~} not to affect the equations of motion. In the most conventional formulation the source terms have been 
anticipated as 

iSw[A-,pA] = Jd''xtT{pA{x)W[A~]{x)}. (3) 

with the Wilson line along the light-cone temporal direction. 



W[A^]{x) ^Vexp 



ig dx+A-{x+,x)T- 



(4) 



defined in the color adjoint representation. Here the source matrix is defined in the adjoint basis; pA = P'aT"- . As 
usual V stands for time (i.e. x^) ordering. The choice of Sw[A~ ,pa\^ however, is not unique and there are as many 
variants as one likes that satisfy gauge invariance and lead to the same equations of motion. Recently it has been 
revealed in |19| that a careful treatment of the fast gluon integration gives rise to a logarithmic form, 

iSw[A-,pa] = J d''xtr{pA{x)\iiW[A-]{x)}, (5) 

which was proposed first in motivated by Wong's equation. We shall adopt the latter |SJ| rather than the 
conventional choice Q throughout this paper when the explicit form is necessary, though this choice should be 
irrelevant to the final results at least within the scope of the JIMWLK problem. 

Now let us see what follows from the generating functional (Q). The classical picture of the color glass condensate 
results from the tree-level approximation of The expectation value of an operator of gauge fields, (0[A])a, is 
evaluated via {{0[A]))pj^ defined as 

{{0[A]))p^= — [ VAS[A+]0[A]e'^™^'^^+'^'^^^''P''^ (6) 
za J 

with the generating functional za[pa] = J^'^A5[A'^] g'^YM+tSw fQj- states with a given pA- Then one can estimate 
(C'[^])a as an average over the pA-distribution and make a classical approximation; 

{0[A])a = [ DpaWaIpa] mA])),, / VpaWa[pa]0[A[pa]], (7) 



using the solution of the Yang-Mills equations of motion A';^{x) as a function of the source. 



6Sym [A] 



SA^ix) 



5Sw[A ,pa] 
_/ 6A^{x) 



= 0. (8) 

A 
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This is in fact a stationary-point approximation and expected to be legitimate when yyA[pA] allows for predominant 
contributions from large pA ~ l/<? and thus large .A^[pa] in the saturation regime. Here it is understood that the 
weight functional is properly normalized; /I?paWa[pa] = 1- Since pA is just an integration variable in 0, the 
A-dependence in the functional form of Wa solely governs the evolution of the quantity (C'[A])a as A changes. This 
evolution is what the JIMWLK equation describes. 

The explicit classical solution to the QCD equations of motion, D^F'^^ — (5''+pA = (with an ansatz = which 
makes the source term simple), is given in great details in literatures 10]. The covariant gauge {df_iA'^ = where 
we follow the convention that covariant gauge quantities are specifically denoted with a tilde) has a one-component 
solution which is as simple as 



where Q;^(af) is a solution of the Poisson equation in the two-dimensional transverse plane, that is, 



(9) 



\y±)plix ,y_L) 



Att 



In- 



\x± - y_ 



(10) 



Here /i^ is an infrared regulator that eventually cancels in the calculation for gauge invariant objects. The cancellation 
of /i^ naturally occurs; l//i gives a scale of long-range force from the source pA. The force extends to infinity causing 
an infrared singularity at /i — > 0. If an object is color singlet as a whole and its finer structure inside p%{x) is not 
resolved from a distance, there should be no color force any longer at distance and thus no infrared singularity. 

The solution in the light-cone gauge (yl+ — 0), which is necessary in the calculation of quantum fluctuations, is 
available by the gauge rotation A^ — > V{A'^ — d'^/ig)V^ eliminating A'^, that results in 



with the gauge rotation matrix 



^^a[PA]=-6^'^^V,-, 



(11) 



exp 



«5 



dz a%{z 



(12) 



Here V represents the ordering in unlike x^-ordering in Q). In this paper, in order to avoid tangled notations 
like pf^, we would use only p\ with a remark about the gauge choice if necessary. In the above discussion p\ is all 
in the covariant gauge, while the light-cone gauge solution (|11|) satisfies the equations of motion for the source in the 
light-cone gauge, namely, T4-,_ooPaKJ- written in terms of the covariant gauge source. 

From these relations we can choose to use Q;^(a;) to express a function of p^(x). Since p%{x) has the support over 
^5 ^ l/^i so does a\{x). One might seem perplexed at this at first glance; since the ai^-structure is such 
localized, a^(x) carries a large longitudinal momentum \p~^\ > A like fast partons generating p/^{x). At the same 
time the classical fields should describe wee partons with < A. This is not a contradiction, however. The point 
is that the scale separation in p+, as we have mentioned, makes sense not in the covariant gauge but in the light-cone 
gauge only. Actually ^'[pa] in the light-cone gauge is delocalized and takes a nonvanishing value at x~ > 1/A. [It 
should be noted that the retarded boundary condition is implicit in the choice of H12|l so that only positive x^ is 
taken.] Recalling a coincidence between the momentum and space-time rapidities, the light-cone gauge solution y^*[pA] 
distributing over x~ > 1/A surely corresponds to slow partons. 

We are now ready to write the JIMWLK equation down in a well-known form in terms of a" (x) using the rapidity 
variable r defined in 0; 



d_ 



S 



6a^{xT , x±) 



where the derivatives act only on the edge of the support x ~ x^ and 

1 fd'^zj_ 



IT J (/TT)"^ 

with the transverse kernel, 



ICix^,y^, z^) 1 ~ V\z^)V{y^) - V\x^)V{z^) + V\x^)V{y^) 



IC{x±,y_^,z±) ^ {27rf{z^\^\xi_) ■ (z^\^\y^) = 

a, a, [xj 



z±)^{yj_ - z±_)^' 



(13) 



(14) 



(15) 
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Here V^x±_) = yj_ - ^+00.-00 with the definition UH). 

The derivation of the JIMWLK equation is our goal pursued in this paper. Usually one derives the JIMWLK 
equation from the fact that the generating functional is independent of the separation scale A as we will explain 
in some details in App. ^ We shall here take a different starting point: The important observation pointed out in 
|19|| is that the generating functional is written equivalently in a different form with an auxiliary variable p as 



.A 

VpAfA[p] / VAS[A+] exp 



iSvM [A] 



d^xp^{x)A-{x) 



with the density of states A/a [p] defined as 



Ma[p]^ JvpAWA[pA]e^p --L Jd^j:tT{pAix)lnW[-iS/6p]{x)} 



S[p]. 



(16) 



(17) 



One can immediately prove that 1)161) returns exactly to after the p-integration. Because (II6II is mathematically 
equivalent with Q), one should be able to derive the same evolution equation from the A- independence of Z as written 
in H16|) . We are attempting to make sure this expectation in what follows. Let us see as an instant check how the 
equivalent classical picture of the color glass condensate arises. This is simply concluded from a claim that p°'{x) 
is static, which corresponds to an approximation that non-commutativity is disregarded as we will argue later (see 
discussions around ^^). There is then no discrimination between p\{x) and p°'{x) and moreover A/a [p] is reduced to 
VVAiPA] and two representations become trivially identical. The quantum evolution is, however, not so trivial. The 
derivation of the JIMWLK equation based on is the main issue of this paper and elucidated in Sec. IIIII 



III. DERIVATION WITH QUANTUM COLOR CHARGE DENSITY 

We will derive the JIMWLK equation (|13(l with our emphasis on the quantum nature of color charge operators. We 
would remark here again that the derivation presented in this section is our central achievement of this paper. First 
in Sec. IIII .^ we will explain about non-commutativity among quantum color charge density from the field theoretical 
point of view, and then we will make proceed in the concrete calculus in Sec. IIII Bl 



A. Quantum color charge density 

Before addressing the computational steps, let us clarify the meaning of quantum color charge density and explain 
why we are interested in that. In the canonical quantization, in gauge theories, the color charge density operators 
obey the commutation relation at equal time, 

[p'^{x+,x),p\x+,y)\ ^-igr'^^nx+,x)5^^\x~yl (18) 

which one can easily confirm by using the canonical quantization conditions. The commutation relation H18|l holds 
for the source of any (fermionic or bosonic) matter field. In the operator formalism of quantization, H18|l is necessary 
for consistency between Gauss' law and gauge invariance, that is, the secondary constraint in the nonabelian gauge 
theory (i.e. Gauss' law) is first-class if we use Dirac's nomenclature for quantization of singular systems M H the 
theory is formulated in the functional integral formalism, (|18|) is expressed by the expectation values in terms of 
classical fields. In the functional integral formalism, by construction, the expectation value is taken in a time-ordered 
way and ((TH|l is reexpressed as 

{0[p]{p%x+,x)p\x+-r^,y)-p\x+,y)p%x+~r^,x)]) = -igr'\0[p] p^{x+ ,x)) S^^\x ~ y) (19) 

for arbitrary functionals 0[p]. 

In the JIMWLK problem relevant in the dense regime where the color charge density is as large as ^ 1/g, we 
can ignore the right-hand side of Ijiyi) . Then, non-commutativity does not play any significant role in the JIMWLK 
evolution equation (apart from generating non-local vertices corresponding to the expansion of the source terms (2^). 
It is a challenging problem to transcend the JIMWLK equation towards the dilute regime where (|18|l or 1)191) must 
be taken into account. This attempt has been in part successful so far; the expansion in terms of the charge density 
derivative has been improved up to the forth order p3 | and now the general Hamiltonian for all orders has been known 
|25l l26| . while the complete quantization procedure has not been quite clear yet. In the dilute regime the quantum 
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evolution has also been studied in the color dipole picture 27] which can be connected to the color glass formalism 

In the operator formalism the Heisenberg equation with l|18(l describes the evolution, given the Hamiltonian is 
known as a function of the charge density operators. Since the color glass formulation, as we have seen in Sec. ^ is 
naturally understood in the form of the functional integral, it is desired to formulate the evolution equation consistent 
with (|19|l within the functional integral description. Recalling that the color charge operators can be regarded as spin 
operators aligned in color group space, one can anticipate a functional integral representation from knowledge on the 
spin system. That is, the action for the spin system is known to consist of the dynamical and geometrical parts. In 
the language of quantum mechanics, the action is the integral of the Lagrangian which is given as L{q) ~ pq — H{p) 
with a variable q and its canonical conjugate p. The dynamical and geometrical parts correspond to the integral of 
H(p) and pq respectively. In fact the geometrical part has an interpretation as the topological phase which is often 
referred to as Berry's phase. Therefore, some attempts to augment the formalism with Berry's phase in color space 
have been proposed in order to handle (|19|l properly psL |25| . 

Recently, in 19], a different way to deal with non-commutativity (|19|l has been proposed in the functional integral 
formulation. The advantage of that proposed method as compared with the representation with Berry's phase is that 
the necessary modification to the conventional formulation is minimal. There is no need to add Berry's phase into the 
formulation. The point is that p° (x) appearing in lfTB|) rather than pl{x) fulfills ^ and can be regarded as quantum 
color charge density. That is simply because the color charge density operator corresponds to iS/SA~{x) acting on 
the source action exp[zS'vy [^~, p]]. If we choose the source action as in then iS/SA^{x) is replaced by 

coming from the eikonal coupling. It is obvious at this point that non-commutativity is lost if the x^-dependence 
of p'^{x) is dropped. In the next subsection we will see how this x^-dependence is eventually abandoned in deriving 
the JIMWLK equation. This means, in other words, that one would be able to surpass the JIMWLK equation if one 
keeps the full x^-dependence, which is beyond the scope of this present paper. 

B. Renormalization group equation 

We shall elaborate the renormalization group equation here from the fact that the density of states A/a [p] evolves 
so as to keep the whole generating functional Z independent of the separation scale A. That is, Z given in l|16() 
can be written at a different separation scale 6A for arbitrary b. In this paper b is assumed to be small as usual, 
i.e. 6 <C 1 and we are working up to the leading order of Q;sln[l/6] ~ agdr and will finally take the limit of 6 ^ 1 
to acquire the differential equation in t. Then one will reach MbAip] expressed in terms of A/a[p] accompanied by 
contributions from the integration over degrees of freedom inside the scale from A to 5A. Such a relation is nothing 
but the renormalization group equation for A/Aip] that flows according to the scale A or the rapidity r. 

Let us briefly summarize in advance each step for the calculation procedure for the sake of guiding what will 
transpire below. 

(Step-1) We will first start with the generating functional Z written at the separation scale A. We will decompose 
the gauge fields into the slow and fast parts according to the longitudinal momentum p'^ , and we will expand the 
Yang-Mills action up to the quadratic order in semi-fast gluons a'^ which we will define there. 

(Step-2) We will perform the Gaussian integration with respect to and will reach an expression of the generating 
functional Z with the gauge field fiuctuation below the scale bA. 

(Step-3) We will equate the counterpart of the density of states with MbA[p] after the Gaussian integration. In the 
limit of 6 ^ 1 we will acquire the differential equation, that is, the renormalization group equation. 

(Step-4) We will rewrite the renormalization group equation in terms of yVA[pA]- 

In the rest of this subsection we shall closely explain each step in order. In App. ^ for comparison, we will sketch 
the conventional derivation of the JIMWLK equation directly for Wa[pa] in accord with (Step-1), (Step-2), and 
(Step-3) above. 

(Step-1) In order to eliminate degrees of freedom within the momentum shell to be integrated, we shall decompose 
the gauge fields into the semi- fast gluons whose longitudinal momenta lie in the strip, bA < \p~^ \ < A, the classical 
solution, yl(^[pbA] of ((TT|l . and the rest SA^ with < bA, that is, 

A^,^a^,+ S^'^Al [ptA] + S^-6A- . (20) 

We assume that the transverse fluctuation around .4^ is small enough to be negligible. Actually one can show that 
the only current (defined as the conjugate to SA'^) correlators to be retained in the leading- log approximation are 
correlators of the charge density (conjugate to SA~) Therefore the specific decomposition (|20|l is validated under 
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the present approximation. It should be noted that from the argument below (|12() we know that the classical solution 
appears as a result of the stationary point approximation in evaluating the generating functional, and thus the classical 
part in the above decomposition H20|) should be that for pf,A (not for pa)- This is because is integrated out and 
y^^[pf,A] should correspond to the stationary point of fields which have < bA. 

Then, it is straightforward to expand the Yang-Mills action in the generating functional at the scale A in terms of 
a^' up to the quadratic order. 



//•bA 
VpAf/,[p] / VSAS[SA+] exp 



tSYM[A[nA] + SA]-t d''xp''{x)SA-{x) 



X / VaSla^] exp 

JbA 



i I d'^x 



5S' 



YM 



SA^ix)SAUy) 



jd^xp^{x)a-{x) 



where the derivatives of the Yang-Mills action are taken at A\ [pbA\ + SA'^ . 

(Step-2) Since (|21|) is quadratic in one can easily perform the Gaussian integration to have 



(21) 



^ Va5[a+] exp -i j d'^x d'^y [a^'^ix) ^ i j d 
X |<(y)-z jd'z'Gll'Ay,^'){^ 



5S' 



YM .A- 



SA^^z) 



5^-p^{z) ) Go\^,(z, x) ) G.'^ix, y) 



lab / 



JYM 



= det 



Go-r(^,2/) 



d'^d''y{^^-^'"p''i^)]GZbi^,y) 



\SAii{x) 

1 



6S- 



YM 



SAUy) 



exp 



d^xd^" 



( 6S- 



6S- 



YM 



where the propagator inverse for the semi-fast gluon is written as 



iG^^^\x,y) 



5A>i{x)5Al{v) 



{ [D.D^r'a,. - [D.D^Y' + 25/-^F;,} 6('\x - y) 



with 



D-=d--ig5A-, D' = d'-igA\ F+, = -D.SA' , F^i=d+Ai, F,j = 



(22) 



(23) 



(24) 



in the light-cone gauge. We note that all quantities listed in p4(l are matrices in the color adjoint basis. 

In the rest of this step let us evaluate the determinant in the above expression (|22|) , namely, the determinant of the 
propagator inverse which is expressed in the (a~, a*) basis explicitly as 



*Go-r(^,y) = 



-d+Df + 2gP^''{d+A''A 



-Dfd+ - 2gf'"'^{d+A1) -[D\D^r5,j - [A-Dj]"'' 
from which one can immediately write the determinant of matrix blocks as 



(25) 



det 



Go-r(^,2/) =deth9+^<5 



+2 rabi 



X det 



[DxD^fH,^ - [A^'j]"'' - {-Dfd+ - 2gr'"^{d+A1)]{~d+'^)-^{-d+Df + 2gf'' ''(9+^1^^ )} 



. (26) 



Here the first part, Aet[~d^'^5°'^] 2, is just an irrelevant number, and after some algebra the second part simplifies 
as 



det 



[D^^D^rS,, - 2grn ^[d+A]) - {d+A!i)^ \ 2gT^f' '\ {d+A^)^A 



'^g^r'f'''\d+Ai)^{d+Ai) 



2 racd rdc'b J ac 



1 



1 



d+ 



(27) 



Although this complicated expression looks impracticable to proceed further, we can drop many terms under the 
leading- log approximation as we will explain now. Let us focus on the second term ^ {di/d^){d^Aj) (the first term 
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in the first curly brackets) for the moment to take an example. If one remembers that the determinant originates 
from the terms sandwiched by the semi-fast gluons a{^, one readily understands that the derivative di/d~^ picks up 
the momenta associated with the semi-fast gluons, which is proportional to with 6A < |p+| < A. The terms 

suppressed by can be neglected because no logarithmic enhancement ^ ln[l/&] arises from those terms. Because 
the first term ^[D\D^]'^''Sij contains a term proportional to p^, the logarithmic enhancing terms could appear 
from the part with no such suppression if the determinant is expanded in the power of Let us next 

turn to the second curly brackets in (|27|) which contains Ai{l/d~^){d^Aj) ^ AiAj apart from the suppressed term 
Ai{d^Aj){l/d^) that is to be dropped by the reason we now explained. This would potentially produce logarithmic 
enhancing terms if expanded. Interestingly enough, however, there is a complete cancellation by the contribution 
from the last term in (|27|l . i.e. {d^Ai){l/d^'^){d^Aj) ^ AiAj apart from the suppressed terms. Consequently we 
only have to retain the first term in H27|) and thus arrive at 



det 



(28) 



within the leading-log approximation. Now we make a further approximation on this. By the argument given in [2a | 
the interaction region in which both pA {x) and 5A^ [x) are nonvanishing is presumably small in the two-dimensional 
space spanned by x"*" and x~ , that means both and p~ are substantially larger than the transverse momentum 
scale. Then we can ignore didi as compared with 28^8" inside the determinant. Using d'^SA'~ ~ which is valid in 
the light-cone gauge in the leading-log approximation {5A~ has \p^\ < bA and so d~^SA~ leads to contributions cx b 
which are neglected) the determinant is thus 



det 



det[2D+i:>~] ^ = det[d+ - iga] ^ det[d- - igS A']' 



(29) 



Here we once moved to the covariant gauge to avoid the complexity associated with entanglement between 9* which is 
neglected and A^ in the light-cone gauge. It should be mentioned here that the determinant of the covariant derivatives 
as in l|29|) is gauge independent as it is individually. Thus we could rotate SA~ in the covariant gauge back to SA^ in 
the light-cone gauge without touching anything on det[(?+ — iga]~^, namely, det[9'^ — igSA~] = det[9+ — igSA^]. In 
this paper we do not discuss but just discard det [9+ — iga] which would not play any role in the JIMWLK problem 
but might be important if we consider the evolution of not only the target but also the projectile. 

The fact that det[9^ — igSA^] is gauge independent implies that this be reduced to an irrelevant number if 6A^ is 
eliminated by an appropriate gauge rotation. In other words, the determinant is expanded as a sum of all the Wilson 
loops which turn out to be irrelevant numbers because they sit only on the one-dimensional space. If we require the 
periodic boundary condition in a;"*", which could be suggested from the physical situation of elastic scattering, we can 
consider the non-trivial Wilson loop wrapping over a;"*" whose starting and ending points are connected by periodicity. 
Therefore we can expect that under the periodic boundary condition det[9~ — igdA~] is given as a function of M^[(5yl^] 
defined in Ql. 

The explicit evaluation of (|29|l is possible in the same way as discussed in leading to the the classical color 
charge density as 



det[d- - igSA-] ^= JvSpawa[Spa] exp --^ J d^xtT{SpA{x)lnW[SA-]{x)} 



(30) 



This is a gauge invariant expression, as we have already noted. Here the correction to the classical color charge density 
Spa{x) distributes in 1/A < < 1/6A because the determinant is taken in the space spanned by the semi-fast gluon 
which has bA < < A by definition. Usually only the positive x~ distribution is chosen in accord to the retarded 
boundary condition which is considered as physically relevant. The combinatorial weight function associated with 
the spatial distribution is denoted by wa[Spa\- This part, however, has no logarithmic enhancement ^ ln[l/6] and is 
to be discarded in the leading-log approximation. It should be mentioned here that the term containing lnT4^[(5j4^] 
turns out to be a shift operator acting on logarithmically enhanced terms and thus we should keep it. 

(Step-3) Now that we have accomplished the Gaussian integration, we can identify the counterpart corresponding 
to A/fcA [p] as a function of Ma [p] , which leads to the renormalization group equation for Ma [p] ; 



VpAfbA[p]e~'^''"^P°^^^'^^<^^^ = / VpAfA[p] exp[AS[A[pbA] + SA-,p]\ 



X / VSpA exp 



I d^xtT{5pA{x)lnW[iS/6p]{x)} 



-i f d^x {x)5A^ {x) 



(31) 
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where AS* is the rest of the Gaussian integration (|22|l given by 

AS[A[p,^] + SA-,p] ^ -1(4^^ S'^-P^'i^)) GZ,ix, y) - 5^-p\y) ) • (32) 



(Step-4) Then the evolution equation for the weight functional Whip a] follows immediately from as 

Vp e-^/'^'=^P°(^)'5^a(^) jvpbK 

- y^A[Pbh\ e-5^/'''^*''{'5'^A(^)l"^[-*'^/'5''l(5)} gAS[^[pi.A]+<5^-,p] ^-^J <f'xt,{-pt,(x)\nW[-i&/Sp](x)} 

Here we defined the integration j Vp^A = J VpA J VSpA, and the weight functional is extended as [2^ 



(33) 

5[p\ - 0. 



WAkA]=>VA[pA] n '^['^/'a]- (34) 

l/A<a:-<l/bA 

We leave 5A^ in AS" as it is (not replaced by —i6/6p) on purpose to avoid confusion in the approximation we will 
make later. 

In the dense regime where pA takes a large value, we can make an important approximation here, i.e. we can expand 
the Wilson line term as 

- ^ Jd''xtr{pAix)lnWH5/Sp]ix)} 

where the higher-order contributions involve more derivatives 6/dp, so that they are suppressed when pA or p is large. 
[p is of the same order as pA as we will see shortly.] Here V stands for time ordering as in The first term 
in 1)35(1 placed on the exponential leads to the shift operator which makes f[p] be f[p — pa]- The second term is 
vanishing under the second-order perturbation that the terms up to the quadratic order in SA^ or S/Sp"^ are kept; 
the derivatives in the second term should act on the function of p'^{x) arising as a result of the operation by the first 
term. Then the delta functional shifted by Pa(^) leads to p"'{x) — p^ix) which is x+ independent. Once /0°(a;) loses 
the a;+ dependence, then only the first term in the expansion 1(35(1 remains, which signifies that non-commutativity 
becomes completely irrelevant. Therefore we insist that the expansion ((35() is the key equation to understand how the 
non-commutativity is disregarded in the JIMWLK problem in the high density regime. In Add. 1X1 we will comment 
on the connection of ((35(1 to an approximation made in the conventional derivation of the JIMWLK equation. 
In evaluating ((33(1 . keeping only the first term in the expansions like 1(35(1. one can simplify it as 



'j:,pp-Kfd*^p''ix)SA^(x) ip- 



'JbA 



S[p-pA-SpA]^0. (36) 



One can readily complete the p-integration here, then replace 5A in AS* by the derivative in phA = Pa + SpA, and 
move the derivative onto yVA[p6A] by the integration by parts. Finally we get the evolution equation in the form of 

WbAiPbA] - e^^[-^['^''Al^.V^P.A.PA]>v^[^^^] ^ 0. (37) 

It should be noted that the pbA-derivative hits all pbA from the far left, that is, it acts on pA in AS* as well as on phA in 
VVAiPbA]- A subtle point is that AS' contains not only pf,A but also pA- This is because the origin of pA in AS is traced 
back to the eikonal coupling ^ p'^cl~ for semi-fast gluons. The integration over leads to the renormalization SpA 
for soft gluons, but not for semi-fast gluons themselves. It might seem then that taking the limit of 6 —> 1 in AS is a 
subtle business. This does not matter in fact, however, because there remains no contribution from the pbA-derivative 
hitting PbA inside AS as we shall confirm explicitly below. [We will see another reasoning for that in App. IXI ] 

The rest of this section would be devoted to illustrating the evaluation of AS[A[pbA] — iS/SpbA, Pa]- At this point 
we realize that AS'[^[pbA] — iS/Spi,A, pa] is the same quantity calculated in Sec. 6 of j2^, while ((32() was different. We 
can borrow the results from ^26] but shall try to make our calculations in a self-contained fashion. We will work in 
the light-cone gauge, which eventually leads to the same formula as the calculation in the covariant gauge in of 
course. By using d~^6A~ — one can easily see 



6S 



YM 



dAil 



S^'-pA = S^''2^g^+[r'^^Al[pbA]} + S^'SpA (38) 



AlptA]~iS/SptA ^ ^PbA 
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in the light-cone gauge. One can also express it in terms of the quantities in the covariant gauge by means of 



r 



be 



iNr 



6A 



-A' 



bA 



9 



bA 



(39) 



Here we substituted (fTT|l and rotated the yOf,A-derivative by T4-._oo to make it in the covariant gauge and used a 
formula Vab = {l/Nc)tr[T'^VT^V^. Also we performed the integration by parts to move 9' onto S/SptA- Although we 
skip the detailed explanations, we can use the free propagator, i.e. the propagator in the absence of the background 
field Al^{x), if we employ the light-cone gauge. This is, roughly speaking, because in the light-cone gauge A^^{x) has 
no singularity but only a discontinuity at a;~ = which does not affect the propagator jl81] . The free propagator can 
be further approximated in the interested region as jl^ 



^oabi^^ y) 



An 



(40) 



where dr = ln[l/6]. The other components, Gq^j,, Gq^'j,, and Gq^"^, are suppressed by l/p"*" as compared with Gq^j, 
and we drop them in the leading-log approximation. Therefore, under this approximation, we can neglect 6^~SpA in 
H88|l because no logarithmic enhancement results from the terms with S'^^SpA in AS*. 

Because the first term in (|38|) is a total derivative in and the approximate propagator H40|l does not depend on 
x~ , the integrations over and are almost trivial, so that only the surface terms at x~ ,y~ — ±cxd remain. For 
instance, the integration of the first term in (|38|l with respect to x~ results in 



Sp^{Xr ,x±) 



6p^{0,xj 



(41) 



In the above we have referred to the separation scale by using the rapidity variable r instead of A or bA and have 
used the fact that Pt{x) has a support < a;^ ^ a^r • The matrix V is V^- _^ ~ V+oo,--oo with the definition H12|l . 
The integrand does not depend on time and thus the x^ and integrations are absorbed in the redefinition as 



^ dx'^ 5 1 5 pr{x) 



27r 



5/5pr{x) (see (6.8) in 2^]). After all AS" takes a form of 



d x± 



((9' 



Sp^{xt,x±) 



Sp^{Xr ,X±) 



(0 



6p''{0,Xj_)J\ Sp''{xr,x±) 



Sp^{Xr ,Xj_ 

Vabix^) + (a' 



(.. 



dpb{0,X^) 



Sp^{0,x^] 



6p'^{0,XA 



(42) 



which is nothing but (6.5) in |26j |. It is apparent from this expression that the pT-derivative does not act on pr in AS" 

because (5Vjj/(5p°(a;-) cx {T'')acVX = and (5Vjft/(5p^(0) cx V^^{T'')cb = 0, etc. Then in exactly the same way as in 
|26|. we can change the variable from p'^{x) to a^(a;) by H10|) and with the help of relations. 



(5<(0,a;i) fe^(a;r,a;_L) 



Vlix^) = Vabix^) 



5a^^{xT ,a;±)' 



(43) 



it is immediate to confirm that we finally reach the JIMWLK equation H13() . Let us focus on the second term in 142(1 
to take an example. Using ((10|l . ((14|l . and ((43|l in order we have 



dr 
2^ 
dr 
2^ 
dr 



d^xAd 



Sp'^iXr ,XA 



VlbM 



= -— I d^xj_d^y^d^z^{z^\-^\x± 



J yt( ) 

Sa-{xr,x^)Sab{0,y^) "^"^ ^' 
J d^x^d^y^d^z^lC{x^,y^,z^) ^ , ^cb{y ±) , , ^_ , Vl{z±), 



l{Xr,Xj_) 



5a''^{xr ,yx) 



(44) 



which is identical with the second term inside the square brackets in the evolution kernel p4|l . 

Before closing this section, we shall comment upon the meaning of (|43|l . We can easily validate these relations 
supposing that it is only V or V'^ on which the ar-derivatives act. It does not mean, however, that YVtIpt] is a 
function of V and V'^ . In fact, the derivatives can be moved from on >V^[p^] to on ©[pr] on averaging over a,- (or pr) 
in O . Thus, it is sufficient to require O [pr] written in terms of V and V'^ , which is actually the case if we consider the 
scattering between the color glass condensate and the projectile consisting of partons in the eikonal approximation. 
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From the practical point of view we can say that the relations of H43() would not lose generality of the evolution 
equation. 

This treatment may well arise a question about the ordering problem. It is obvious from 144|l that 5Vcb{v±)/ Sa'^ ix±) 
is nonvanishing and so 5/Sa'^{x~) and S/Sa'^{0) do not commute after H43I) is applied, while they should be commuta- 
tive by construction. This is not a paradox, however. For example, in H44() . there is certainly a non-zero contribution 
from SVcbiy ±)/Sa'^{x±), but such a term is precisely compensated for by the term arising from non-commutativity 
among ofr-derivatives acting on V or , that is, 



be 



1 



(45) 



if these derivatives are supposed to act on V or Therefore, there is no ordering problem at all unless the quantities 
in l|43() are split apart; in H43(l an ordering like {S /Sa'^{y ^)){S /Sa'^{x±))Vcb{y ±) is not allowed, which generates extra 
terms and causes the ordering ambiguity. In other words, the ordering realizing in the JIMWLK equation (|13(l is 
determined uniquely in a way that H45|l would bring about no ordering problem. 

In the present calculation we did not treat the virtual term cr° (see (|A9|)) in an explicit manner. The virtual term 
is, in fact, included implicitly in the correct ordering as pointed out in 26] and it will be clearer in the discussions in 
App.EI 



IV. CONCLUSIONS 



We have derived the JIMWLK equation in the functional integral formalism. Our starting point is different from 
the conventional derivation; we used a different representation of the QCD generating functional consisting of the 
density of states and the simple eikonal coupling. In this way, we could clarify how the non-commutative nature 
of quantum color charge density becomes irrelevant in the approximation that the color charge density is large. 
Actually the expansion of the Wilson line in the non-local source term has turned out to be responsible for loss of 
non-commutativity. We further discuss in App. some connections between the approximations made in the new 
derivation and their counterparts in the conventional derivation. The comparison between two different methods 
would give us a deeper understanding of those approximations. 

The next step would be how we can improve the approximation to handle non-commutativity in the functional 
integral formalism. The commutation relation between the functional derivatives with respect to the classical gauge 
field (see (|45f) ') originates from that the derivatives act on the gauge invariant Wilson lines, and there is no difficulty 
in this case. In contrast, in the case of the commutation relation between the color charge density operators, the 
functional integral description for that is non-trivial because the color charge density in the functional integral is not 
an operator like the derivative but a given classical field. From our analysis we would presume that the non-local 
source term written in terms of the Wilson line takes care of this, that is, if one is capable to take full account of the 
source term, one could reach a general evolution equation in terms of classical quantities alone. In our formalism such 
a manipulation would be possible, at least in principle, by diagonalizing the classical color source, and rotating it back 
in color space. We are making progress in this direction. The point is, once it will be accomplished, the quantization 
is automatic and there is no need to impose the commutation relation between color charge density like in |26|. 

Although it is far from obvious whether the renormalization group equation is closed in a color glass condensate 
formalism in an intermediate region between the dense and dilute regimes, it would be a challenging problem to attack 
it in the functional integral formalism. We believe that the present formulation would provide a useful information 
for that purpose. 
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APPENDIX A: COMPARISON TO THE CONVENTIONAL DERIVATION 



In this paper we have elaborated the derivation of the JIMWLK equation in the framework with the source terms 
consisting of the density of states and the simple eikonal coupling. Now, let us look over the conventional derivation 
to grasp some connections between two different strategies. 
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The derivation we address here is not faithful to the original one in literatures. We shall rearrange the calculations 
and put our emphasis on some aspects which are related to our derivation in the text, though the main building blocks 
are not quite new. The weight functional >Va[pa] changes according to A in a way to keep the generating functional 
Z intact, from which the evolution equation follows. That is, Z can be written at arbitrary separation scale 6A as in 
(P). We shall take (Stcp-1), (Step-2), and (Step-3) respectively as in the derivation we addressed previously. 



(Step-1) The gauge fields are decomposed into the semi-fast gluon with &A < < A, the classical solution 
^a[PA + ^p\i and the rest with |p^| < 6A , and then the action, S = Sym + Sw, is expanded up to the quadratic 
order in terms of a{j'; 



VAS[A+] exp iS[A, 



bA 



X exp 



VdAd[SA+] / VaS[a+] exp iS[A + SA,pA 

JbA L 



d^2 



5A'i{x) 



SAil{x)5A^,{y) 



(Al) 



It should be mentioned that in this case the action contains Sw [A~ , pa] given by (|3Jl that generates non-local terms 
in a;+. This is an important difference from what we have seen in our derivation in Sec. IIIII 

In the dense regime where SA^^ is small (|Alll can be further expanded in terms of (5A^ to read the quantum 
corrections to the source as 



j\?aS[a+]exp -1 j d^x d^y ai^ix) G-l'^\x,y) al{y) - i J d^x 6p-{x){-a- (x) + SA' (x)) 



where we have used the equations of motion including the source terms, i.e. 

SS[A,pa] 



SA^{x) 

and we have defined the corrections to the source, 
Sp'^ix) = Sp^^^'^ix) + (5p(2)"(a;) = - [ d^y 



-A[pA+Sp] 



Sp^ix) 



(A2) 



(A3) 



SAa{x)SA-^{y) 



<iy) 



d yd z 



5^S 



SAa{x)SAUy)SAi{z) 



al{y)a^^{z). 

(A4) 



This approximation that one keeps only the linear term, —i5p"'{x)5Ag^ (x), corresponds to the expansion of (|35|l and 
is responsible for loss of non-commutativity. In writing (|A2p we have defined the propagator for with Sw included 



(A5) 



(that makes a difference from Gq^^ appearing in Sec. 1111 as 



_PS_ 

'x)SA 

where the argument pA of G~l°''^ refers to the source contained in Sw- 



iG-l'^'>[pA]{x,y) = 



5A':,{x)5Al{y) 



(Step-2) The generating functional with the separation scale A is now rewritten in a form from which we can 
identify the weight functional at scale &A, that is, 



VpaWa[pa] / I?M(5[M+]e'^[-^+*^'^^l 



VpaWa[pa] I VaS[a'\e ? 

bA 



+ 1 ld''xd''ya'^{x)G-i'"'ix,y}a"^{y}+iJd-^xSp''(x)a-{x) 



bA 



VSAS[5A+]e 



iS[A+SA,pA+Sp] 



, (A6) 



where the last term in ljA2p gives the renormalization of the source pa ^ PA + 5p in 5* in the above expression. We 
shall introduce an auxiliary variable p^A which is constrained as ~ Pa ~^ ^P- It should be worth mentioning that 
(5/9 is a function of a'^ and should be interpreted as the quantum color charge density in a sense in discussions in 
Sec. nil Al In the dense approximation, however, as we have explained before, 5A~ is taken only up to the first order 
and at that stage we have already lost the discrimination between quantum and classical color charge density. By 
introducing phA and integrating over pA, we reach 

r-A 



(A7) 



VpbA VpaWa[pa] / I?a(5[a+]e^5-f'''"'''«<("'^M.'"'(-'2')<to)+^/^'-*p"(-)''^(-)(5[pfeA-pA-M 

J JbA 
JbA 
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where pA inside G^l°'^[x, y) must be altered according to pA PbA — Sp, though this could be ignored in the Gaussian 
approximation with respect to a^{x). Otherwise, the already quadratic term —^a^G~^°'''a'j^ could produce higher 
order terms. 

Here we shall remark that the role of the 5p°'{x)a~ (x) term should deserve further investigations. If we treat it in 
the Gaussian approximation, it cancels Gz]^'' and G~^^ , and is irrelevant to the final results. It has been anticipated 
in [2^, however, that the complete integration over a'^{x) including the full Yang-Mills action would result in the 
Wess-Zumino term for 5p°'{x). In this paper, we shall just drop 5p°-{x)a~{x) hereafter. This treatment is valid in 
deriving the JIMWLK equation in the dense regime. 

Expanding the above expression in terms of 5p (with the 5p°'{x)a^{x) term dropped), we have 

>V.A[P-.A] = >VA[m] -/.^^X^(>VAbA]."(x)) + \ J 

where the correlation functions are defined as 



a'^ix) = {Sp^'^'^ix)) = ld*yd*z ..... . -Gl^iy^z) (A9) 



1 (d^yd^z ^ 

2 J ^ 8Aa[x)bAl[y)8Al{z) 

X^\x,y) = (<5p(i)'^(:r)<5p(i)^(y)) = / d^d^zd^w ■ ^. /'f.,. , • G^,^(z,^) (AlO) 

J 8Aa{x)SA>^{z) (y)M2(u)) 

within the Gaussian approximation. This is exactly the same expression as (3.53) in 16] and the evaluation of these 
correlation functions goes just in the same manner. We shall comment that in this case the determinant as a result of 
the Gaussian integration would not lead to the determinant 1)30(1 because the propagator defined in this section does 
not contain 8A~ in it. 

As we have already mentioned, we could have ignored the difference by bp inside the propagator within the Gaussian 
approximation. Then the pf,A-derivatives in l(A8p would not hit the pA dependence in the propagator. This fact seems 
to correspond to the arguments in Sec. IIII Bl that the pbA-derivatives vanish acting on pf,A in AS'. We cannot see the 
mathematical correspondence transparently, however, for the propagator l|A5p involves non-local contributions from 
Sw{A~ ,p\. 

It is known that the explicit evaluation reveals an interesting relation, 

(All) 

from which ((A8|) simplifies as 



2i «»(!/) 



y^AM .ilMyjl^,->(..y)^^WM. (A12) 



This corresponds to the fact, as we mentioned before in Sec. IIII Bl the correct ordering implies the virtual term cr°(a;). 
We obtained in Sec. IIII 51 the same ordering as in l|A12p which is determined by the condition that the commutation 
relation (|45|l becomes free from the ordering problem. The ordering as shown in HA12(I implicitly signifies the presence 
of the virtual contribution through (|Alip . In this way, we should consider that our derivation contains the virtual 
term as well in an implicit way, which has also been mentioned in part in p^. 
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